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ABSTRACT 

We develop a method to infer or rule out the presence of an atmosphere on a tidally-locked hot super 
Earth. The question of atmosphere retention is a fundamental one, especially for planets orbiting M stars 
due to the star's long-duration active phase and corresponding potential for stellar-induced planetary 
atmospheric escape and erosion. Tidally-locked planets with no atmosphere are expected to show a 
Lambertian-like thermal phase curve, causing the combined light of the planet-star system to vary with 
planet orbital phase. 

We report Spitzer 8/im IRAC observations of GJ 876 taken over 32 continuous hours and reaching 
a relative photometric precision of 3.9 x 10~ 4 per point for 25.6 s time sampling. This translates to 
a 3-cr limit of 5.13 x 1CU 5 on a planet thermal phase curve amplitude. Despite the almost photon- 
noise-limited data, we are unable to conclusively infer the presence of an atmosphere or rule one out on 
the non-transiting short-period super Earth GJ 876d. The limiting factor in our observations was the 
miniscule, monotonic photometric variation of the slightly active host M star, because the partial sine- 
wave due to the planet has a component in common with the stellar linear trend. The proposed method 
is nevertheless very promising for transiting hot super Earths with the James Webb Space Telescope and 
is critical for establishing observational constraints for atmospheric escape. 

Nutzman fc Charbonneaul (120081) 



1. INTRODUCTION 

Super Earths are a recently discovered class of exoplan- 
ets, loosely defined to be 10 M® or less. Since these plan- 
ets have low masses, they likely consist substantially of 
rocky material, making them the first analogs of terrestrial 
planets in our solar system. The search for super Earths 
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siting super Earths ( Rouan et a l. 2009) anticipated in the 
next few years. 

A fundamental question about rocky planets is whether 
or not they have an atmosphere. The atmosphere can reg- 
ulate the temperature and protect the surface from harm- 
ful solar or cosmic rays, issues central to surface habit- 
ability. Extreme ultraviolet (E UV) radiation, wind s, and 
coronal mass ejections (CMEs; iLammer" et alJl2007ft from 
the host star are the main mechanisms for driving atmo- 
spheric escape. M stars have a much longer active phase 
compared to solar- like stars, during which a star's EUV, 
wind s, and CMEs are strong (jScalo et al.ll2007tlWest et all 
2008). The issue of atmospheric retention, therefore, is es- 
pecially relevant for super Earths in sho rt-period orbits 
around M stars (e.g.. ILammer et al.ll2007h . 

We are motivated to observationally determine whether 
or not a tidally-locked close-in super Earth orbiting an M 
star has an atmosphere. The goal is to observe the planet 
and star in combined light and search for thermal phase 
variations peaking at the substellar point, indicative of a 
bare rock planet with no atmosphere. 

The idea to study transiting hot supe r Eart hs in 
this way was briefly mentioned in SelsiJ (|2004f ) and 



same idea for non-transiting terrestrial planets orbiting at 
the Earth's semi-major axis from a sun-like star. The idea 
would be to use a Ter restrial Plane t Find er /Darwin- type 
interf erometer (e.g., ICockell et al.l 120091 : lLawson et al.l 
2008, and references therein) to null out the star light 
to directly image the pla net over the course of its orbit. 
iGaidos fc Williami (|2004f ) explored light curves of planets 
with and without atmospheres and oceans, for a variety 
of obliquities, eccentricies, and viewing geometries. Com- 
bined light phase curves of about half a dozen hot Jupiter 
transiting planets h ave been obtained with the Spitzer 
Space Telescope (e.g v lHarrington et alj[2006t ICowan et alj 
l2007t iKnutson et atll2007l 12009ft . Even though our first 
proposal to observe a hot super Earth phase curve was in 
2005, three attempts at observations were needed to reach 
a reasonable S/N and to mitigate instrument systematics. 

We begin in §2 with a description of the motivation and 
background to detect an atmosphere on a tidally-locked 
super-Earth. In §3 we describe Spitzer IRAC observations 
of the planet-hosting M dwarf star GJ 876A and report 
on upper limits to inferring an atmosphere on the short- 
period super Earth G J 876d. In §4 we discuss prospects for 
the method to infer or rule out an atmosphere on tidally- 
locked super-Earths, including transiting exoplanets. 

2. MOTIVATION AND BACKGROUND 

The motivating factor for observations to discriminate 
between a bare rock and a planet with an atmosphere is 
simply that basic atmospheric escape estimates and calcu- 
lations are inconclusive because they are highly dependent 
on assumed planetary paramets (such as mass, primordial 
atmosphere mass, and the past EUV radiation history of 
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the host star). 

Atmospheric escape is a complex process thought to 
be driven by nonthermal escape processes or by EUV ra- 
diative heating. The high temperatures and consequent 
expanding upper atmosphere enables thermal escape of 
gases. To illustrate the main uncertainties in EUV heat- 
ing, we use the energy-limited escape framework (Lecave- 
lier des Etangs 2007). The gravitational potential energy 
of the planet [J] is 



E p = — 



GM p m atm 



(1) 



where G is the gravitational constant, M p is the planet 
mass, R p is the planet radius, and j3R p is the radial extent 
of the planet exosphere. Here m atm is the mass of the 
planet atmosphere, which we assume to be low enough 
to ignore the radial atmosphere structure in our estimate. 
The incident EUV power [J s _1 ] incident on the planet is 

P uv = nR 2 p F EUV . (2) 
We assume here that the radius where the EUV energy 
is absorbed is close enough to the planetary radius. Tak- 
ing the ratio of equations ([1} and ^ , and an active-phase 
average EUV flux < Few > leads to the lifetime of the 
initial planet atmosphere 



G 



R p 



(3) 



As long as the temperature is high enough, the absorbed 
stellar radiation is approximately instantaneously reradi- 
ated on a bare rock planet. We can see this in a first order 
way by comparing the radiative flux to the conductive flux. 
The radiative flux is, 



aT 4 , 



(4) 

where a is the radiation constant a — 5.670 x 
10~ 8 J K~ 4 m~ 2 s _1 . Conduction is in general very inef- 
ficient in rocky material. Heat conduction is defined by 



Fn = 



dQ _ k dT 
dtdA dz ' 



(5) 



0e < Few > 

Energy- limited escape assumes that some fraction, e of the 
incident UV flux heats the atmosphere and drives escape. 
Note that (3 > 1 and e < 1. This equation is valid provided 
that the upper atmosphere is heated strongly enough for 
the relevant elements (e.g., H, or C) to escape. 

We can see the main uncertainties in atmospheric ther- 
mal escape from equation ([3]), even though the equation 
is approximate at best. First, if a planet is transiting we 
know the planetary density. If the planet is not transit- 
ing and we assume the planet is solid, the density varies 
by less than a factor of five for reasonable interior com- 
positions (Mercury-like to a "dirty water" planet) an d 
for a range of planetary masses (e.g.. ISeager et al.ll2007f) . 
More uncertain is the initial outgassed planetary atmo- 
sphere mass, which could vary widely, from less than 1 
percent of the planet's total mass up through 5 percent 
and eve n as high as 20 percent for an initially water-rich 
planet ( Elkins-Ta nton fc Seagerll2008| ). An equally large 
unknown is the incident EUV flux. All stars have very 
high EUV flux during their so-called "saturation phase", 
where levels can reach thousands of times present-day so- 
lar EUV levels. The high EUV is most critical for planets 
orbiting M stars: M star saturation phases can last 1 to 
3 billion years for early M stars and 6 to 8 billion year s 
for late M stars (M5 to M8 respectively; IWest et al1l2008f) . 
If the initial planetary atmosphere has not survived dur- 
ing the star's saturation phase, a new atmosphere could 
develop in the star's quiet phase if the outgassing rate is 
higher than the escape rate. 

The basic idea for observations is that a tidally-locked 
planet with no atmosphere, i.e. a bare rock, will show a 
thermal phase variation that can be observed in the com- 
bined light of the planet-star system. We assume short- 
period sup er Earths are tidally locke d based on dynamical 
timescales ()Goldreich fc Soteilll966|) . 

Q is heat in the form of kinetic energy of the motion of atoms and molecules. 



where Q is healQ, A is the cross-sectional heating area, 
dQ/dtdA is the heat flux, k in units W m _1 K _1 is the 
thermal conductivity. For Earth's surface materials, k ~ 3. 
Conduction describes how heat is transfered from a hotter 
to a colder region, as indicated by the negative tempera- 
ture gradient in the above equation. 

We can estimate the conductive flux by adopting differ- 
ent temperature gradients dT j dz. A tidally-locked planet 
has a maximum temperature gradient across the surface 
created by the drop off of incident stellar flux away from 
the substellar point. This drop off follows (1 — cos^ 1 / 4 ), 
where 6 = is defined at the substellar point and 9 is 
the angle away from the substellar point surface normal, 
and using equation ((4]). The temperature change across 
10 m, 100 m, or even 1 km is much less than a fraction 
of a percent. The same line of reasoning can be applied 
to regions away from the substellar point, also resulting 
in the dominance of radiation over conduction. Hence the 
conductive flux across the surface is essentially negligible, 
and radiation back to space will dominate over conduction 
away from the heated surface element. 

We can also estimate the conductive flux from the sur- 
face down into the planet interior. We can conserva- 
tively take Earth's subsurface temperature; a more mas- 
sive planet is likely to be hotter in the subsurface due to 
more heat-generating radioactive decay, and tidal heating 
may also play a role. The Earth's subsurface temperature 
at a depth of 500 m is about 293 K. This is deep enough 
so that the primary heat source is radioactive decay of ele- 
ments in the Earth's interior. If we take GJ 876's substellar 
temperature of ^650 K, then, using equations ([4]) and ([5]), 
the radiative flux is 5000 times greater than the conductive 
flux, showing that radiation back to space will dominate 
over conduction of energy down into the planetary inte- 
rior. Even if we take an artificially extreme case of a 10 K 
m temperature gradient just beneath the planet's sur- 
face, the radiative flux is still a few hundred times greater 
than the conductive flux. We note that Earth is hotter 
in the interior than at the surface and has a temperature 
gradient of about 0.02 to 0.03 K m" 1 ; on Earth there is 
no conduction of absorbed sunlight down into the crust 
because conduction only operates from hotter to colder 
regions. 

In summary, as long as the temperature is high enough, 
reradiation will dominate over conduction. We may say 
that if a hot super Earth is hot enough for its atmosphere 
to have escaped, it should also be in a regime where rera- 
diatio n dominates over conduction. iGaidos fc Williams! 
( 2004) have used a different line of reasoning to argue that 
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the thermal inertia of even cooler rocky planets is negligi- 
ble. 

For instantaneous absorption and reradiation, the plan- 
etary thermal phase curve for a tidally-locked bare rock 
takes the sam e expression as that of a Lambert sphere 
(lSobolevlll975h . 

&a — — [sin a + (ir — a) cos a] . (6) 

The phase angle, a, is the star-planet-observer angle. 
With this definition, a = 0° corresponds to "full phase", 
a > 170° corresponds to a thin crescent phase, and the 
planet is not at all illuminated at a = 180°. The phase 
angle in terms of inclination is 

cos a = — sin i sin 2-7T0, (7) 

where <p is the orbital phase. For the phase-dependent 
thermal infrared flux ratio of a bare rock we hav e an ex- 
pression similar to that of re flected light (e.g., ISobolevI 
ll975t[Charbonneau et al.l ll999) with the geometric albedo 
A g replaced by (1 — A g ), 

F P {<t>,i) n , v (R P \ 2 1,. v , 

— = (1 — Ag) I — - I — [sin a + {it — a) cos a\ . 

■l ^ V Cb J 7T 

(8) 

We are therefore looking for a sinusoidal variation like 
a Lambert sphere, depending on the inclination of the 
planet. Detection of a Lambertian thermal phase curve is a 
necessary but not sufficient condition to identify a tidally- 
locked exoplanet without an atmosphere. Some tidally- 
locked exoplanet atmospheres, including those with thin 
atmospheres or with strong absorbers at high altitudes 
in thick atmospheres will shows a thermal phase curve 
that peaks on the planetary day side. The only way to 
discriminate between a tidally-locked bare rock exoplanet 
and tidally-locked planet with a thin or thick atmosphere 
is that atmospheric winds likely shift the planetary hot 
spot away from the substellar point. The data needs to 
be good enough to pinpoint the hottest point on the plan- 
etary day side to see whether or not it is consistent with 
the substellar point. 

3. OBSERVATIONS AND DATA ANALYSIS 

We observed GJ 876d, because at the time of our 
first proposed observations (2005) GJ 876d was the 
only short-period super Earth candidate available for the 
method to infer th e presence of an atmosphere. GJ876 
([Rivera et alj2005h is a dM4 star with T ef f = 3350±300 K 
(jManess et al.ll2007ft . At a distance of 4.69 pc, GJ 876 is a 
nearby star and very bright in the near-IR (K = 5.0). 

The planet GJ876d has Miami = 5.89 ± 0.54M©, 
P = 1.93776±7x 1 0~ 5 , and a semi-major axis of 0.02 AU 
(jRivera et al.l l2005h . A non-transiting planet, GJ876d's 
orbital inclination is unknown. D espite being a rela - 
tively slow rotator (96.7-day period; I Rivera "eTatl 12005). 
and relatively inactive, GJ876 has a low level brightness 
variability of order 0.05 magnitude at visible wavelengths 
(|Shankland et alJl2006ft . 

We observed GJ 876 at 8/im wavelength using 
Spitzer/IRAC. The observations were nearly continuous 
for 32 hours, starting on July 16, 2008 at 01:50 UTC, us- 
ing an IER in two segments. We acquired a total of 4204 
data cubes in IRAC subarray mode, each data cube com- 
prising 64 exposures of 0.4-sec duration, and 32x32 pixel 



spatial extent. IRAC observations at 8/j.m are subject to 
an increasing sensitivity as the detector is exposed to light, 
an effect termed the "ramp" (Doming et al.l l2006f ) . We 
used a "pre-flash" strategy to flood the detector with high 
flux levels and thus force the ramp to its asymptotic value 
prior to observing GJ876: we observed the compact HII 
region WC89 for 59 minutes immediately prior to observ- 
ing GJ 876. 

Our analysis used the Basic Calibrated Data produced 
by version SI 8. 0.2 of the Spitzer analysis pipeline. We ap- 
plied the corrections recommended by the Spitzer Science 
Center (SSC) to correct for variations in pixel solid an- 
gle and flat-fielding for point sources. We corrected the 
effect of energetic particle hits by applying a 5-point me- 
dian filter to the 64- frame time history of each pixel within 
each data cube. We set the threshold of the median fil- 
ter to equal 6% of the peak value in the average stellar 
image, and we corrected discrepant pixels to the median 
value. We performed aperture photometry on each of the 
64 frames within each data cube, using a circular aperture 
with a radius of 4.25-pixels. We centered the aperture for 
each frame using a parabolic fit to the brightest 3 pixels in 
the X- and Y-profiles of the star, each profile derived by 
summing over the orthogonal coordinate. We varied the 
aperture radius to minimize the scatter in the final photo- 
metric time series. We measured a background value using 
an annulus with inner and outer radii of 6 and 12 pixels. 
Each photometric point (Figure 1) was constructed by av- 
eraging the frames withi n each data cube, om itting the 
1st and 58th frame (see lHarrington et al.l 120071 ) . We also 
omitted frames that were internally discrepant by more 
than 4(7 from the other frames within each data cube. 
We subtracted an average background value for the entire 
time series (not frame-by-frame) because this minimized 
the scatter in the time series photometry. 

Figure 1 (top panel) shows the photometric time series. 
We found that a very small ramp effect (about 0.1%, not 
illustrated on Figure 1) was still present for the first ~ 2 
hours of the data. We simply dropped the first 159 min- 
utes of data from the analysis, keeping the last 3828 data 
cubes (orbit phases > 0.3). 

The most prominent feature on Figure 1 is a 2-millimag 
stellar flare that occurs near p hase 0.355. Flares are quite 
common on active M-dwarfs (|Reid fc Hawlevl 120051 ) , al- 
though none have been previously reported on GJ876. 
However, three facts indicate that this is a real flare: 1) 
Spitzer photometry for solar- type stars has a proven stabil- 
ity (e.g.. lKnutson et al.ll2007f ). with no instrumental effects 
of this type, 2) the event has the characteristic rapid rise 
and exponential decline characteristic of M-dwarf flares, 
and 3) Knutson et al. (2007, and private communication) 
observed a similar event on the M-dwarf companion to 
HD 189733. We fit and remove a flare model from the pho- 
tometry. The flare model has a linear rise and exponential 
decline, requiring 4 parameters (rise rate, time of peak, 
amplitude, and decay rate) and we fit for the minimum y 2 
using a gradient-expansion algorithm. We emphasize that 
the flare is not a significant impediment to our analysis, 
because it has a well defined shape and occupies only a 
relatively small fraction of the light curve. 

In addition to the flare, we observe a gradual decline 
in brightness that is equivalent to 0.6 millimag per day. 



4 



GJ876 is known to be variable (|Weisl [l994'). and ground- 
based photometry in the Johnson B band by P. Sada 
(private communication) for times bracketing our Spitzer 
observations shows a decrease in brightness at a rate of 
2.0 ± 0.2 millimag per day. This is part of a 97-day pe- 
riodic ity due to stellar rotation, derived by iRivera et alj 
( 2005), and confirmed by Sada's photometry concurrent 
with our observations. After removing the flare, we fit and 
remove the corresponding section of a 97-day sine wave, 
with amplitude scaled down to fit our time series data. 
Note that we expect more muted changes due to stellar 
rotation at this long IR wavelength as compared to visi- 
ble wavelengths. Based on temperatures for starspots in 
active late-type stars (Strassmeier et al. f992), Sada's 2.0 
millimag per day at Johnson B scales to 0.8 millimag per 
day at 8/im, in good agreement with our observed slope. 
This gradual decline in stellar brightness is a more seri- 
ous problem than the stellar flare for putting limits on the 
planet signal, as explained in the next section. 

The middle panel of Figure f shows the data after re- 
moval of the stellar and linear baseline in units of the aver- 
age stellar brightness. The scatter per point is3.9xI0~ 4 , 
only slightly larger than the photon noise (3.54 x 10~ 4 ). 
We thus obtain about 90% of the photon-limited signal-to- 
noise ratio. The distribution of photometric values (Figure 
f, middle panel) after removal of stellar effects is indis- 
tinguishable from Gaussian. Under our hypothesis of a 
planet lacking an atmosphere, the flux modulation due to 
the planet should peak at phase 0.5, computed from the 
ephemeris values given by IRivera et alj (|2005l ). and up- 
dated by Rivera and Laughlin (2008, private communica- 
tion). However, we point out that removal of the baseline 
slope will affect the extraction of a planet signal, since the 
partial sine-wave due to the planet has a component in 
common with a linear baseline. To account for this effect, 
we fit and remove a straight line from the planet model, 
and we fit the difference curve (a "reduced sine") to the 
data in order to estimate the best-fit planet amplitude and 
error. 

A least-squares fit of a reduced sine wave peaking at 
phase 0.5 gives an amplitude of f .26 ± f .71 x I0~ 5 , i.e., 
zero within the errors. To determine the error, we fit sim- 
ilar reduced sine waves to f 00,000 bootstrap Monte-Carlo 
trials based on permutations of the Figure f (middle panel) 
data. Tabulating the distribution of sine-wave amplitudes 
from these trials, we find an excellent fit to a Gaussian 
distribution, with a 3er limit = 5.13 x I0~ 5 . From the 
photometric calibration of our data, we find the flux from 
the star to be 864 mJy, so our 3er limit on the amplitude 
due to the planet is 44/xJy. At a 99% confidence level, the 
flux modulation limit due to the planet is 34/jJy. To illus- 
trate the 3 — a limit in a "chi-by-eye" fashion, the lower 
panel shows a sine wave with 44/xJy amplitude in compar- 
ison to our data binned into 8 phase intervals, as well as 
the nominal best-fit amplitude (not a detection). 

4. RESULTS 

Our Spitzer 32-hour dataset shows no evidence for the 
sinusoidal variation peaking at phase 0.5. To interpret the 
observed limit on flux modulation, we need to account for 
the orbit inclination of the planet and the albedo. We as- 
sume a tidally-locked rotation with no heat redistribution, 



because we are testing the no-atmosphere hypothesis. At 
each assumed inclination value we compute the thermal 
radiance curve of the planet, emitting as a Lambertian 
sphere. Given an assumed inclination, the Doppler data 
yield a mass, and we compute a radius for various compo- 
sitions using the m ass-radius relationships calculated by 
ISeager et alj (|2007t ). This information is sufficient to de- 
fine the amplitude of the planet's thermal emission curve 
at all values of orbit inclination and albedo (see equa- 
tion ©). 

A given upper limit on the planet thermal phase curve 
corresponds to a locus of values in the inclination-albedo 
space illustrated in Figure 2. Tidally-locked atmosphere- 
less planets with albedos and inclinations in the regions 
below the curves in Figure 2 would have high enough ther- 
mal phase curve amplitudes to be detected, and are there- 
fore ruled out by our observations. In other words, if, for 
a given composition GJ 876d has an albedo and inclina- 
tion such that it lies below an upper limit curve in Fig- 
ure 2, we could infer the presence of an atmosphere (imply- 
ing efficient heat redistribution between the tidally-locked 
planet's day and night side). We emphasize again that the 
observations actually rule out combinations of inclination, 
albedo, and radius. By knowing GJ 876d's mass at a given 
inclination, we can translate a mass and adopted interior 
composition to a radius, using planet interior models. 

We are able to make robust (99.8% confidence level) 
upper limit statements that are inclination-, albedo-, and 
composition-dependent, corresponding to regions below 
the curves in Figure 2. As an example, we would infer 
the presence of an atmosphere if GJ 876d is a pure water 
ice planet that has an albedo less than 0.35 if the planet's 
orbit al inclination is 50 degrees (the inferred inclination 
from [Rivera et alj 120051 ). We would infer the presence of 
an atmosphere for a planet with the same composition 
and an albedo less than 0.79 and a n inclination of 84 de - 
grees (the inferred inclination from iBenedict et al.ll2002j ). 
A pure water ice planet is not a particularly useful example 
because planets are expected to have some rocky compo- 
nent. We can also infer the presence of an atmosphere if 
GJ 876d is less dense than a planet with Earth's bulk com- 
position (approximately 32.5% iron core and 67.5% silicate 
perovskite mantle) and has a geometric albedo less than 
0.25 and an inclination of 84 degrees. Because a bare rocky 
planet likely has an albedo lower than 0.25, we would infer 
the presence of an atmosphere on GJ 876d for this compo- 
sition and inclination. For other albedo/inclination values 
for which we could infer the presence of an atmosphere 
for a given hypothetical composition of GJ 876d, see the 
curves in in Figure 2. We are unable to make conclu- 
sive statements about planets more dense than those with 
Earth's composition, because size decreases with increas- 
ing density, and smaller sized planets would have a smaller 
thermal phase variation, rendering them undetectable in 
our data. 

5. DISCUSSION 

We reported nearly continuous 32 hours of Spitzer 
IRAC data of GJ876d. The goal was to infer whether 
or not GJ 876d has an atmosphere, because a tidally- 
locked planet with no atmosphere is expected to show 
a Lambertian-like thermal phase curve. We are unable 
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to either conclusively infer the presence of an atmosphere 
or rule one out. We can present inclination- and albedo- 
dependent statements of what kind of planet atmospheres 
are ruled out if GJ876d is a planet of relatively low average 
density for a solid planet (see Figure 2). 

The limiting factor in our data analysis is variation of 
the stellar flux. The presence of a slow linear drift in- 
creases the size of an atmosphereless planet that can hide 
in the data. Even though GJ 876d is a quiet M star, it 
still s hows variability o n a rotational time scale of 96.7 
days (jRivera et al.ll2005t ). 

The serendipitous detection of a flare (Figure 1) is fur- 
ther evidence for activity on a very slightly active star. We 
emphasize that the flare itself should not impede efforts to 
extract the planet signal, because the flare has a specific 
shape and a relatively short duration. The detection of 
infrared microflares on M stars may provide a new way to 
study M star activity not available from the ground. 

Ground-based visible photometric monitoring of the star 
during the IR space-based measurements, combined with 
the rotation rate of the star and a star spot model (area 
and temperature of spots) would help to subtract any stel- 
lar variability in the light curve. Additionally, space-based 
IR monitoring before and after the planet observations 
would extend the data baseline to help precisely define the 
stellar variability. Future observations should also use con- 
tinuous observations of the planet-star system, not just for 
maximizing S/N, but also to avoid the detrimental "ramp" 
effects that would be caused by observations of other tar- 
gets of random brightnesses. Monitoring of a planet over 
two or more orbital periods would also help disentangle a 
linear stellar variability component. 

Higher S/N would also mitigate effects of stellar vari- 
ability. Our upper limit was degraded because of overlap 
between a component of the planet's variation and the 
stellar variability. We fit to the data with a "degraded 
template" . Having more photons would help to get more 
precision, and tighter limits than using such a degraded 
template. We caution that with more photons we might 
encounter another aspect of stellar variability at a lower 
level. Our observations may be the first to indicate the 
difficulty even low levels of stellar variability have on ob- 
servational studies in the combined-light of the planet and 
star, for any kind of exoplanets orbiting an M star. 

Knowing the inclin ation of GJ 876d w ould help in in- 
terpreting our data. iRivera et alj (|2005l ) provide an in- 
clination for the GJ 876 3-planet system of ~ 50 degrees, 
based on stability arguments for the two ^Jupiter-mass 
planets GJ 876b and GJ 876c, assuming that the lower 
mass GJ 876d is also coplanar. Hubble Space Telescop e 
Fine Guidance Sensor observations ([Benedict et alj |2002) , 
in contrast to the theoretical simulations, give an incli- 
nation of GJ 876b of 84±6 degrees. Even with improved 
observations and theoretical simulations, the inclination of 
GJ 876d may remain unknown, unless we adopt the dy- 
namical assumption that it should be co-planar with its 
more massive planet siblings. 

Transiting planets provide a much better opportunity 
for observations to infer whether or not an atmosphere ex- 
ists from a thermal phase curve compared to a planet with 
unknown orbital inclination. Because a transiting planet's 
orbital inclination is near 90 degrees, the planet-star flux 



contrast is maximized for a bare rock planet. The radius 
of a transiting planet is known, enabling a specification 
of the albedo without resort to mass-radius interior mod- 
els. (While non-transiting planets with known inclination 
would still help with interpreting the results, the unknown 
radius remains a complicating factor.). We can use Fig- 
ure 2 to understand what kind of limits one could put for 
comparable data for a transiting planet. If GJ 876d were a 
transiting planet at 90 degrees inclination with no evident 
thermal phase variation, we would infer an atmosphere on 
GJ 876d for almost any interior composition. The excep- 
tion is planets denser than a planet with a Mercury-like in- 
terior composition (approximately 70% iron core and 30% 
silicate mantle); slightly higher S/N data would be needed. 
For a planet with a Mercury-like composition, the albedo 
upper limit is 0.15. Most bare rocky bodies in our so- 
lar system have albedos lower than this value Mercury's 
geometric albedo is 0.1 and the Moon's is 0.12). One ex- 
ception is Io, Jupiter's atmospher eless moon, with a high 
geometric albedo (0.63 ICoxll2000l ) due to fresh solid sul- 
fur deposits from episodic volcanic outbursts, themselves 
caused by tidal friction with Jupiter. Although many of 
the icy satellites of the outer plan ets also have very high 
geometric albedos (0.63 lCoxll2000l ). an ice-covered surface 
is not expected for hot planets like GJ 876d and others 
at very close planet-star separations. Data with S/N ~3 
times higher than our data set reported here would rule 
out such high albedos even for the smallest plausible plan- 
ets. 

We have discussed how to infer the presence of an atmo- 
sphere on an exoplanet. Can an atmosphere be ruled out 
based on detection of a thermal phase curve? In principle 
a planet with a thi ck atmosphere having inefficient heat 
redistribution fe.g.. lHarrington et alj|200a ) could produce 
a signal similar to a bare rock. Measuring the phase of 
the orbit variation to high precision is needed; a bare rock 
must peak very close to phase 0.5, but a planet with a 
thick atmosphere should have a detectable phase shift of 
the thermal maximum, even for relatively inefficient heat 
transport. Discriminating between a planet with a very 
thin atmosphere and a bare rock may be more difficult, 
and requires further study. 

Observation of the combined light of the planet and star 
as a function of phase may also reveal the presence of a rel- 
atively strong planetary magnetic field. For short-period 
planets, the planet and star magnetic fields could interact 
to produce a phase curve peaking at phases near 180 de- 
grees with a period equal to that of the pl anet's orbit, as 
has been detected for a few hot Jupiters (IShkolnik et al.l 
120051: iDonati et alJl2008tlShkolnik et al.ll2008t) . The aues- 
tion of super Earth magnetic fields is a serious one. With- 
out a protective magnetic field the atmosphere might be 
further eroded by stellar wind and high energy particles. 
A detection of a relatively strong magnetic field would in- 
dicate the presence of a partially liquid interior, aiding 
interpretation of interior composition (and even further 
interesting for short-period planets under strong tidal in- 
teraction.) 

We note that Mercury is not tidally-locked to the sun, 
(i.e., in a 1:1 resonance). Instead Mercury, is in a 3:2 res- 
onance, rotating three times for every two orbits about 
the sun. Models show a low probability for a planet on 



6 



a circular orbit to become trapped in the 3:2 spin-orbit 
resonance compared to the 1:1 resonance, unless it goes 
through a chao tic evolution of an ecce ntric orbit as Mer- 
cury likely has ([Correia fc L askar 2004). Planets on circu- 
lar orbits are therefore more certainly tidally-locked than 
eccentric planets which may have had the opportunity to 
be captured into higher order spin-orbit resonances. 

The James Webb Space Telescope, with its 58 times 
more collecting area than Spitzer will be able to make high 
S /N observations of short-period super Earths. M stars are 
good targets, because they are bright at IR wavelengths 
and the planet-star radius ratio is high. 

In conclusion, inferring an atmosphere on a short- 
period, tidally-locked super Earth that lacks a thermal 
phase curve is possible for a transiting planet of any in- 
terior composition and albedo, provided that the data set 
has higher S/N by about a factor of 3 than ours. During 
the JWST era we anticipate the birth of studies of atmo- 
sphereless bodies and the concomitant understanding of 



atmospheric escape. 
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Fig. 1. — Upper panel: photometry of GJ876 for 3828 data cubes, each representing 62 images. The initial small detector ramp, represented 
by 376 data points (159 minutes) was omitted from the analysis and is not shown. The red line is a fit of a linear (0.002 per day) brightness 
decrease of the star, plus a flare of amplitude 0.002 near phase 0.355. Middle panel: Data from the upper panel after the fit (red line) is 
subtracted. Lower panel: Residuals from the middle panel binned into 8 phase intervals, and shown in comparison to our 3cr upper limit on 
the flux modulation (solid curve, amplitude = 5.13 X 10 — 5 , 44 ^tjy) peaking at phase 0.5 as expected for a planet without an atmosphere to 
redistribute heat. The dotted curve is the nominal best fit (not a detection). 
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Fig. 2. — Inclinations and albedos for which an atmosphere can be inferred. Each line is the locus of orbit inclination and albedo values 
that would produce an 8 fim sinusoidal flux modulation corresponding to our 99.8% and 99.0% confidence upper limits (44 & 34^tjy). Color 
encodes planet composition: pure water ice planets (solid blue curve); "dirty" ocean planets (dashed blue curve, 6.5% by mass iron core, 
48.5% silicate mantle, and 45% ice outer layers, or different combinations with the same mass and radius); an Earth-like composition (red 
curve, 32.5 % iron core and 67.5% silicate mantle); and a Mercury-like composition (green curve, 70% iron core and 30% silicate mantle). A 
GJ 876 with no atmosphere is allowed by our data above each line; the regions below the line are forbidden at that confidence level. In other 
words, if GJ876d has a composition, orbit inclination, and albedo lying in the regions below the curves, then wc would infer the presence of 
an atmosphere. 



